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c I ! 1. Introduction 

^"'i The typical theoretical prediction for an extragalactic neutrino flux at E v ~ 10 16 — 10 18 eV 

is the generic upper bound on the diffuse neutrino spectrum from AGN jets obtained by 
Mannheim et a In their model AGN source is not completely thin for a CR flux, and the 
normalization of the resulting CR spectrum is such that the CR intensity does not exceed 
the proton spectrum obtained from observations. At E v ~ 10 16 eV this bound is 
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A review of the recent achievements in high energy neutrino physics and, partly, neutrino astrophysics 
is presented. It is argued that experiments with high energy neutrinos of natural origin can be used for 
a search of new physics effects beyond the electroweak scale. 



j{E)E 2 < 3 x 10~ 7 GeV • cm" 



(j(E) is the differential diffuse neutrino spectrum). More conservative bounc corresponds 
to the case when the sources are completely transparent for CRs: j(E)E 2 < 2 x 10~ 8 
(independently on the neutrino energy if j(E)ocE~ 2 ). 

At higher energies, E v ~ 10 20 eV, extragalactic neutrino flux can have several compo- 
nents: the guaranteed component is the "cosmogenic" neutrino flux resulting from interac- 
tions of CRs with relic photons of CMB(cosmic microwave background). According to the 
most optimistic prediction^ the value j(E)E 2 for the cosmogenic flux may be as large as 

3 x 10~ 7 GcV • cm" 2 • sr" 1 • s" 1 

at E v ~ 10 20 eV. Other components of neutrino flux at ultrahigh energies are hypothetical: 
e.g., neutrino from decays of topological defects, neutrino from interactions of new (exotic) 
hadrons with mass ~ (2 — 5) GeV in CR ftuj0with relic photons etc. 

Experimental bounds on the extragalactic diffuse neutrino flux are on the level 

j{E)E 2 re 10~ 6 GeV • cm" 2 • sr" 1 • s" 1 

obtained, at E pa 10 15 eV, in experiments on neutrino telescopes, and, at higher energies, on 
large air shower arrays (see, e.g., review 



i 



2,2008 18:20 WSPC/Guidelines-IJMPA Bugaev26oct 



2 E. V. Bugaev 

2. Neutrino-nucleon interactions 
2.1. uN DIS in Standard Model 

Generically, the differential cross section for the process viN — > IX at large neutrino 
energies strongly depends on the behavior of parton distribution functions (pdfs) at small 
values of Bjorken variable x and large values of Q 2 . Besides, the gauge boson propagator 
effects are sufficient (leading to a nonlinear rise of <j v n{E u )) when Q becomes comparable 
to the electroweak scale, Q 2 ^S> M^. In the energy region which we are interested in, 
E v > 10 7 GeV, one has 

Qchar « Mw, ^char < 10" 4 - 10" 5 . 

In QCD-improved parton model perturbative QCD corrections to the nonperturbative 
pdfs (taken from experiment) are calculated using two approaches: DGLAP scheme (re- 
summations of the perturbative expansion retaining leading terms in ln(-nr — )) and BFKL 

QCD 

scheme (resummations retaining leading terms in ln(i)). Both these logarithms become 
large in the relevant kinematic region. Since, very roughly, x-dependence of pdfs of sea 
quarks has a power-law behavior, oc x~ 03 , for x -C 1, it follows that 

a vN {E v ) cx E n v 3 . 

The power-law rise of the cross section implies a violation of unitarity at very high ener- 
gies. 

It is rather evident that there must be non-linear ("higher twist") QCD corrections to 
the standard calculation, those which are beyond the linear evolution DGLAP and BFKL 
schemes. Physically, these non-linear corrections arise due to the growth of parton (gluon) 
density in the nucleon target at small x. Taking into account the recombination gluon ef- 
fects leads to a taming of the fast rise of the cross section, in agreement with unitarity 
requirements. 

There are different approaches accounting these gluon screening effects. To feel the 
order of magnitude of corresponding corrections to a V N{E„) it is enough to compare the 
results of two models: unified BFKL/DGLAP + screening model^ and phenomenological 
colour-dipole modelPl n the latter model, the DIS is considered in the laboratory system: 
virtual weak gauge boson at some distance from the target fluctuates into qq-pair inter- 
acting with the target via multiple gluon exchanges. This approach, in its original form, 
does not take into account QCD evolution at all. It appears that at extremely high neutrino 
energies, E v w 10 12 GeV, the difference in the predictions of these models for the total 
cross section does not exceed factor 2 — 3 (see Ref. 6 ). 

The ct u n{Ev) scales, approximately, as E® 363 for E v > 10 16 eV. Using this one 
can easily obtain the estimate for the rate of contained events in neutrino telescopes.^ 
According to this estimate, the value j(E)E 2 should be larger than 

10" 7 GeV • cm" 2 • sr" 1 • s" 1 

to be measurable at E v > 10 19 eV at 1 km 3 detector, if the SM prediction for the total 
cross section is correct. 
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2.2. TeV scale gravity models for a 

It is easy to show that the gravitational scattering of any two particles (via graviton ex- 
change) becomes strong if y/s ^> Mp\ w 10 19 GeV, i.e., at inaccessibly high energies of 
colliding particles. 

In TeV-scale quantum gravity models proposed by Arkani-Hamed et there are two 
most important independent parameters: Mb (fundamental gravitational, or Planck, scale) 
and n - the number of extra (compact) dimensions. The value Md is connected with the 
4-dimensional Newton constant through the compactification radius R. If the compactifi- 



cation volume, V = 2itR, is large, Md, which is equal to I -y^ ) , is much smaller 



than Mpi. In TeV-scale gravity models one has Md ~ 1 TeV. 

In these models, because of the compactification, the extra n components of the gravi- 
ton momentum are quantized. To an observer in the usual 4-dimensional space-time, the 
graviton would appear to be a massive Kaluza-Klein (KK) state. The result of the summa- 
tion over KK states in the expression for the amplitude of gravitational scattering is 



where A(w Md) is the cut-off constant (the Born amplitude is divergent but after eikon- 
alization the amplitude becomes finite). We see from this expression that the amplitude is 
large if s/s > Mo sa ITeV, rather than y/s > Mpi, as in the 4-dimensional theory. Due 
to s 2 -dependence, the gravitational scattering amplitude becomes, at some energy, larger 
than any other amplitude. 

The gravitational amplitude and cross section can be calculated reliably if the process 
is semiclassical, i.e., when the Schwarzschild radius (Rs) of the colliding system (which is 
the classic value) is much larger than quantum lengths, X p (it is higher dimensional Planck 
length) and As (de Broglie wave length of the colliding particles). It is easy to see, that 
for the semiclassical picture of the gravitational scattering the inequality y/s S> Md must 
take place (the "trans-Planckian energy regime"). In this regime gravitational interaction 
dominates over other gauge interactions and is a semiclassical process. 

The important characteristic of the gravitational scattering is the impact parameter b. If 
b ^> Rs one has — - <C 1 (t is the square of the momentum transfer) and it corresponds 
to a small scattering angle limit (and elastic scattering). If, in opposite case, b < Rs, the 
collision leads to a production of a black hole with Schwarzschild radius Rs and mass 
Mbh equal, roughly, to cms energy yfs. 

The condition yfs ^S> Md corresponds to the existence of the minimum value of BH 
mass which can be reliably produced in the collision: 



a is the model parameter. The cross section of the BH production in vN collision is equal 





BH — 



min 



V s min = aM D , a > 1, 
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where fi(x, fi) is pdf for the parton of type i, /i ss Rg , and a^H \ s me cr0 ss section of 
BH production in i^-parton collision at cms energy y/xs, 

a™(xs) m kRUMbk = y/xi)Q{y/xi - aM D ). 

The differential elastic gravitational scattering of neutrino on nucleon, j^, is calculated 
(using multiple graviton exchangJIII) as a function of y = E »~ E » , This cross section 
grows as y decreases. The small y region corresponds to long distance processes where 
neutrino interacts with a parton and transfers only a small portion of its energy, surviving 
after interaction. 

2.3. Neutrino interactions in string theory 

If cms energies of collisions are close to the fundamental scale of gravity, i.e., if </s « Md 
("Planckian region"), the classical description cannot be trusted. String theory provides the 
best hope for understanding the regime of strong quantum gravity, and for computing am- 
plitudes at energies close to Md- So, the models with extra dimensions must be embedded 
in realistic string models, in which the unification of gravity with the SM takes place. 

In string theory the massless graviton is the zero mode of a closed string, whereas the 
gauge bosons of SM are the lightest modes of an open string. The scattering amplitudes 
in string theory are amplitudes of string exchanges and are described by the formulas of 
Veneziano's (exchange of an open string) or Virasoro's (exchange of a closed string) type. 
Correspondingly, one expects the presence of string Regge (SR) excitations of the graviton 
and gauge bosons (analogously to resonances in Veneziano model). By duality arguments, 
if there are amplitudes with SR excitations exchanges in the t-channel, having correspond- 
ing t-channel poles, these amplitudes will also have s-channel poles. 

In lepton-quark scattering the s-channel resonances are leptoquarks. From the reso- 
nance amplitude 

v + q -> X;l -> v + q 

(J is a spin of the resonance) one can obtain the cross sectionp]cr^(i/q) = a (yq — > X£), 
and in the narrow width approximation one has 

u J n {vq) ~ 9 H (s - nM 2 s ) , 

where g is the gauge coupling constant and Ms is the string scale. One can assume that 
Ms ~ Md ~ 1 TeV, and, in this case, the SR resonances are "TeV strings". So, we have, 
in fact, the GUT unification at TeV scale. 

2.4. Search for a new physics 

Neutrino experiments with neutrino of natural origin can be effectively used for a search 
of a physics beyond the standard model. There are different possibilities. 
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Study of contained events in large neutrino telescopes. At large enough energies of 
neutrino, say, E v > 500 — 1000 TeV, most part of such events may be due to gravi- 
tational elastic scattering of neutrino on the nucleon target or BH production by neu- 
trinos. One can study the angular dependence, up-down ratio as well as the energetic 



distribution 



12 



2. Study of through-going muons with E > Eq, produced in BH production processes in 
large neutrino telescopesU^ 

3. Study of quasi-horizontal air showers which occur at rates exceeding the predictions of 
SM and have distinct characteristics (see, e.g., Ref. 14). Although greatly reduced by 
hypothetical BH production process (the cross section of which can be, in TeV gravity 
models, about 10 4 nb at 10 20 eV), neutrino interaction lengths, L — 1.7 x 10 7 km w.e. 
(^), are still far larger than the Earth's atmospheric depth, which is L = 0.36 km w.e., 
when traversed horizontally. Neutrinos therefore produce black holes uniformly at all 
atmospheric depths. As a result, the most promising signal of BH production by CRs 
is quasi-horizontal showers initiated by neutrinos deep in atmosphere. Distinct charac- 
teristics of the showers are following: anomalous electromagnetic componentP^large 
multiplicity; large ^ ratio; they look nucleus-like (for not a small X max ); they have a 
curved front, with particles well spread in time. 

4. Earth-skimming idea. The large detectors of air showers can register the fluxes of UHE 
neutrinos and, simultaneously, will be able to measure <7„at at energies as high afl^ 
10 20 eV. Several proposed experiments plan to detect UHE neutrinos by observation 
of nearly horizontal air showers (HAS) resulting in atmosphere from v-air interactions. 
At E « 10 20 eV the a vN in SM is about 10~ 31 cm 2 . The air shower probability of 
HAS production is proportional to cj u n- In addition to HAS, the experiments can also 
observe up-going air showers (UAS) initiated by muon and tau leptons produced by 
neutrinos interacting just below the surface of the Earth. The expected rate of UAS de- 
pend non-linearly on the cross section (if o v n larger than 10~ 32 cm 2 , UAS probability 
is inversely proportional to ovtv). Thus, by comparing the HAS and UAS rates, the 
cross section can be determined. Taken together, these rates ensure a total event rate, 
weakly depending of the value of <j v n- If there is a new physics (e.g., BH production), 
this idea is also useful. Really, a large rate of quasi-horizontal showers may be at- 
tributed to either an enhanced neutrino flux or an enhanced BH cross section. While an 
enhanced flux increases these rates, a large BH cross section will suppress them, since 
the hadronic decay products of BH evaporation will not escape the Earth's crustP^l 

3. Neutrinos and UHECR puzzle 
3.1. Z -burst model 



It had been proposed^ that the primaries which propagate across distances above the 
GZK zone (w 50 Mpc) are neutrinos, which then annihilate with relic neutrinos within this 
zone to create a flux of nucleons and photons with energies above Eqzk- The annihilation 
cross section is relatively large (ps 10~ 5 mb) near Z-resonance and the resonance neutrino 
energy depends on the neutrino mass. 
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There are three main difficulties of the Z-burst model. 

1. Primary protons have to be accelerated to extremely high energies, E > 10 23 eV, in 
order to produce in astrophysical sources (via pp and pj reactions) UHE neutrinos. The 
photons produced in the same reactions have to be absorbed inside the source otherwise 
the diffuse background of MeV-GeV photons will be too large. Hidden sources are 
unable to produce neutrinos of such high energies because the decay length of pions 
becomes equal to its scattering length. The luminosities of sources required in this 
model are too highpU (10 45 - 10 47 ) erg/s. 

2. Even if one assumes that neutrinos originate from decays of superheavy dark matter 
particles through the only channel, X — vv, the problems with diffuse photon back- 
ground are not avoided because higher order corrections to this process give rise to 
electroweak cascades transferring large fraction of energy to photons and electrons. ^ 

3. The required UHE neutrino fluxes in the Z -burst model are very large, they are almost 
excluded by the new experimental limits from FORTE and GLUE experiments (as well 
as by the new limit from EGRElEO). 

It was shown in the work of Gelmini et a/EHI using the AGASA data, that the nonob- 
servation of CR events at E CR > 2 x 10 20 eV implies a lower bound w 0.3 eV on the 



neutrino mass. Since this value exceeds y Ato 2 ^ from neutrino oscillation experiments, 
the bound applies to all three neutrino masses. If there is neutrino mass hierarchy then it 
follows from SK data that the mass of the heavier neutrino is ps 0.04 eV. It is argued by 
Gelmini et at, that AGASA data are incompatible with such a low value because in this 
case the resonance energy is too high and this predicts too many CRs beyond the AGASA 
end point, 2 x 10 20 eV. So, the bound leaves only a small interval for neutrino mass around 
« 0.3 eV (if Z-burst model is correct). 

The flux requirement obtained in the work of Gelmini etai, is about 7x 10 -36 (eV-m 2 - 
sr • s) -1 , and such a high value is marginally excluded by GLUE and FORTE experiments. 
The relic neutrino density could be enhanced by either gravitational clustering or by a 
large lepton asymmetry. Calculations argue against of such a clustering of relic neutrinos 
if neutrino mass is so small (m„ w 0.3 eV): the overdensity S of neutrinos in our Local 
Group of galaxies is < 10, on a length scales > IMpc (see Ref. 21). Also, BBN physics 
gives the bound on a lepton asymmetry, < 0.1, from which it follows that the extra 
contribution to p v from degeneracy is very small |22l 

3.2. Neutrinos as CR primaries 

It was noted above that within the SM the neutrino-nucleon interaction cross section at 
10 20 eV is about 10~ 4 mb, i.e., on five orders of magnitude smaller than necessary to 
produce air showers starting high in the atmosphere. It appears also that models with extra 
dimensions predict cross sections which on a factor « 10 2 larger than that in SM. It is 
clearly not enough for explanation of the GZK puzzle. 

It was argued recently ^3 that the cross section at E > 10 18 eV can be greatly (on a fac- 
tor 10 5 — 10 6 ) enhanced by nonperturbative electroweak instanton contributions. The pro- 
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cesses induced by electroweak instantons (which represent tunneling transitions between 
topologically inequivalent vacua) violate baryon + lepton (B+L) number and are character- 
ized by the large multiplicity of final state particles (quarks, gauge and Higgs bosons). The 
transition rate of the tunneling process is exponentially suppressed at low energies when 

£cms « 4£ sph « ^— w 30 TeV. 

aw 

Here, -E S ph ~ 8 TeV is the sphaleron energy. At high temperatures and high energies 
the transition rate can be unsuppressed, leading to observable effects in CR experiments. 
According to Ref.|23| at w 10 20 eV the cross section is about 3 mb. 

According to the calculation of Bezrukov et al. based on a generalized semi-classical 
approach, a severe exponential suppression of the cross section is extended up to energy 
30i? S ph ~ 250 TeV. Bezrukov et al. did not calculate the cross section, but estimated only 
the upper limit of the exponential function. Therefore one cannot obtain the reliable cross 
section from their paper although it is tempting to do thisl^H 



4. Neutrino from dark matter 

The existence of dark matter (DM) is the experimental evidence for new physics beyond 
the SM. The most recent WMAP datJ2Hgive the amount of cold dark matter (CDM) as 

0.095 < ficDM^ 2 < 0.129, (2cr C.L.) 

and the total matter density 

0.126 < n m h 2 < 0.143. 

There are many particle dark matter candidates for a stable WIMPs (weakly interacting 
massive particles) which could have the relevant relic density. 

Indirect DM searches have been proposed^ to observe the products of DM annihila- 
tion including neutrinos. WIMPs, which scatter elastically in the Sun or Earth, may become 
gravitationally bound and, over the age of the solar system, they may accumulate in these 
objects, greatly enhancing their annihilation rate. Neutrinos can escape the Sun or Earth 
and can be registered in a detector. 



4.1. SUSY DM 

All superpartners are charged under a discrete symmetry called i?-parity (+1 for SM par- 
ticles and —1 for their superpartners). i?-parity conservation guarantees that the lightest 
superpartner, being odd under i?-parity, is absolutely stable and becomes a DM candidate. 

The desired WIMP is the lightest neutralino x? which is a mixture of the superpartners 
of the hypercharge gauge boson (6°), the neutral SU(2)w g au g e boson (w ) and the two 
neutral Higgs bosons, 

X? = ai6° + a 2 w° + a 3 h° u + aji° d . 
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The recent progress in particle physics and cosmology favors the region of the SUSY 
theory parameters where the neutralino is a Bino-Higgsino mixture rather than Bino-like. 

The rate of neutrino production in WIMP annihilations is highly model dependent. 
Neutralinos are Majorana fermions, and the Pauli exclusion principle suppresses the direct 
production in the annihilations of pairs of light fermions, so only indirect channels of 
neutrino production are available, through the decay of particles produced in the processes 

XX^T + T -,bb, W + W~,ZZ, HH, tt. 

If neutralino is Bino-like, the production of SU(2)w g au g e bosons is suppressed. 

If neutralino exists and is heavy enough and if it is really a mixture Bino-Higgsino, the 
neutrino signal from the Sun direction can be detected by large neutrino telescopes (see, 

fool 

e.g., Ref. ^Sl). The limits on neutralino mass from experiments and from calculations of 
relic number and mass density are very weakEH 

20 GeV <m x < 600 GeV. 



4.2. Kaluza-Klein DM 

There is one class of models with extra dimensions, in which all of the fields of the SM 
propagate in these dimensions, not only gravitons. Each SM field has an infinite tower 
of KK partners with identical spins and couplings and masses of order n/R. Naturally, 
in such models the conservation of the momentum along the extra dimensions and, as a 
consequence, the conservation of KK number, takes place. Radiative corrections break KK 
number down to a discrete conserved quantity, KK-parity: (—1)™, where n is the number 
of the KK level. KK-parity ensures that the lightest KK-partner (LKP) at level one, being 
odd under KK-parity, is stable, and can be a DM candidate. 

In the concrete models the LKP is a neutral WIMP which is a linear combination of 
the first KK mode B\ of the hypercharge gauge boson and the first KK mode of the 
neutral SU(2)w gauge boson. 

The calculation of the relic density of KK DM based on WMAP's results for ficDM 
predicts the value of LKP massl^J 

tolkp ~ (600 - 1200) GeV. 

In contrast with LSP, LKP is a vector particle, and the helicity suppression is absent, 
so, LKP can directly annihilate into neutrinos. Besides, the annihilation into pairs of top- 
quarks is sufficient due to a large value of LKP mass. Top-quarks almost always decay in 
the channel t — > Wb. The W, in turn, decays with equal branching ratios into v e e, v^fi, 
v t t. It is important that both decays, t's and W's (and also r's) take place without any 
energy loss. So,in this model, the Sun is effective source of high energy (E v w 100 GeV- 
1 TeV) tau-neutrinos and (due to oscillations) muon neutrinos. 



4.3. Superheavy DM 

It is well known that the assumption that DM is a thermal relic is too restrictive. Unitarity 
bound on the annihilation cross section and the assumption of thermal equilibrium in the 
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early universe lead to the bouncN^I 

nh 2 > 0.1 (M/10 5 GcV) 2 , 

where f2 is the ratio of the DM mass density to p c today. However, if WIMPs are not 
thermal relics, their masses are not thermodynamically constrained, i.e., they may be even 
superheavy ("wimpzillas"). Moreover, their present day abundance does not depend on 
whether they have strong ("simpzillas"), weak, electromagnetic, or only gravitational in- 
teractions. 

Such particles can be, e.g., gravitationally produced near the end of inflation, due to 
non-adiabatic change of the scale factor during the transition from the de-Sitter to the 
radiation dominated phase. ^ 

High energy neutrinos are produced by the simpzilla annihilations, which produce a 
quark or gluon pair which then fragment into hadronic jets. Neutrino fluxes from the Sun 
were calculated by Albuquerque et a/PH. It was shown that v T and i/ M fluxes can be rather 
large and measurable by neutrino telescopes, with E v in the interval (50 GeV - 1 TeV), the 
spectrum is very similar with the case of KK DM (the channel t — > Wb also works). One 
should note that the typical E v value in the spectrum is much lower than the simpzilla mass 
(M « 10 s — 10 16 GeV), due to the large numbers of neutrinos produced per annihilation. 

Unfortunately,direct dark matter search experiments exclude JJ the most natural values 
of simpzilla mass (those which are comparable with the inflaton mass in chaotic inflation 
models (« 10 12 GeV)). 
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